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Abstract 



We study in a model-independent way anomalous CP-violating tbW effective 
couplings that might arise from new physics in the processes pp — > tW~X and 
CN ! PP — > iW + X, followed by semileptonic decay of t and i. These processes have a 

dependence on effective tbW couplings both in the production process as well as in 
the decay of the t or t. We propose several CP- violating asymmetries constructed 
out of variables in the two processes, including and t polarization, and energy 
and azimuthal angles of the decay particles. We find that it is feasible to probe 
a certain CP-violating combination of anomalous couplings at the per cent level 
at the LHC for centre-of-mass energy 14 TeV and an integrated luminosity of 10 
fb" 1 . 
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1. Introduction: The top quark, with a mass close to the scale of electroweak symmetry 
breaking (EWSB), is generally believed to hold a key to the understanding of the mechanism 
of symmetry breaking and responsible for masses of all standard model (SM) particles. 
While the so-called Higgs mechanism is the dominant scenario for EWSB, it has not yet 
! been established experimentally. The currently operational Large Hadron Collider (LHC), 
which is proceeding full-steam with the search for the Higgs boson, can produce top quark- 
antiquark pairs in great profusion. It will shed light on the details of the properties of the 
top quark, and hopefully, of EWSB. Due to its large mass, the top quark decays before 
hadronization effects, thus preserving polarization information in the decay products. The 
study of the top-quark properties through its production and decay at the LHC would thus 
be of immense significance in arriving at a detailed understanding of the one of the major 
mysteries of nature. 

In addition to EWSB, another important phenomenon which lacks full understanding 
is CP violation (CPV). Apart from having been seen experimentally in mixing and decay 
of K and B mesons, CPV is also essential in understanding the observed baryon asym- 
metry of the universe. In the SM, the only source of CPV is the phase associated with 
the Cabibbo-Kobayashi-Maskawa (CKM) inter-generational quark mixing matrix. However, 
the SM cannot adequately explain the baryon asymmetry pQ. Extensions of the SM are 
needed to understand baryogenesis, and examples of extensions which fit the bill are the 
two Higgs doublet models (THDM) and minimal supersymmetric standard model (MSSM), 
which have additional particles and newer mechanisms of CPV. In fact, in these models, it 
is the CP- violating scalar couplings of the top quark which drive baryogenesis [2]. Thus the 
study of CPV, particularly in the top-quark sector, could shed light on primordial processes 
responsible for baryogenesis. 

Apart from ti pair production, which occurs dominantly through strong interaction, 
single-top production, which necessarily proceeds via weak interaction, also has a large cross 
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section at the LHC [3]-[TT], and has already been seen [12]. Single-top production will 
therefore be able to give us information about the size and nature of the weak-interaction 
coupling of the top quark to the b quark and the weak- interaction gauge boson W, i.e., the 
tbW coupling, as well as the tb (33) element of the CKM matrix. Thus, for example, top 
polarization, which is negligible in the SM in top pair production because strong interactions 
conserve parity, would be large in single-top production because of the chiral nature of weak 
interactions (9]. 

In this work we investigate the possibility of probing putative CP- violating tbW couplings 
in a model-independent way using single-top production in association with a W boson. 
There are a number of proposals for the study of CPV in top-pair production at lepton 
|15j . photon [16] and hadron colliders [17]. However, CPV in single-top production has 
received little attention [18j [19]. The reason can be seen to be two- fold. For one, unlike 
in ti production, where the final state is self-conjugate, in single-top production, a process 
involving t has to be related to another involving t, which is not straightforward. The 
other reason is the low event rate for single-top production expected at Tevatron. Single- 
top production at the LHC can be substantial, and it would be worthwhile attempting to 
extract information on CPV, even though it needs more elaborate analysis than in the case 
of pair production. Recent earlier work on CPV in single-top production at the LHC has 
been in tH~ associated production [18] or in the context of flavor- changing top couplings 
[19], but not so far in tW~ production. This, therefore, is the first time that the possibility 
of studying CP-violating tbW couplings at the LHC is being explored. 

There are three distinct mechanisms for single-top production, viz., a) the t-channel 
process bq — > tq', b) the s-channel process qq' —> tb and c) the tW associated production 
process bg —> tW~ |4J. Although the process a) has the largest cross section, the other 
processes also occur at a significant level. Since each process has a distinctive final state, it 
would be possible to study each one of them separately. 

The process c), in which we are interested, is difficult to isolate due to backgrounds, 
the dominant one being top-pair production. The backgrounds have been studied in [H [5] , 
and kinematic cuts which can clearly isolate the tW channel have been suggested there. A 
similar analysis including NLO effects has been carried out in (SlITJ- Recently, with integrated 
luminosities of 0.7 fb _1 and 2.1 fb _1 respectively, ATLAS JT3] and CMS [H] have presented 
results on a search for the tW signal, where both W's decay leptonically. They have been 
able to put an upper bound on the cross section. 

In a recent work |20j , we examined in the context of the LHC the capability of the process 
c) above, viz., single-top production in association with a W, for providing information 
on possible anomalous tbW couplings in a model-independent way. The processes a) and 
b), while sensitive to tbW couplings, could also get contribution from other processes like 
scalar exchange, and do not permit a model- independent approach. In [20], anomalous tbW 
couplings were found to change the production mechanism, giving rise to top polarization 
different from that predicted by SM. In addition, top decay distributions can also get modified 
by anomalous couplings. Thus, a study of decay distributions in the laboratory (lab) frame 
were shown to carry signatures of anomalous tbW couplings firstly through the degree of 
polarization, and secondly through the contribution to the decay process itself. 

Tevatron provides the only existing direct limits on anomalous tbW couplings through 
W polarization measurements. These results [21] lead to a limit of 0.3 at the 95% confidence 
level (CL) on the magnitude of the tensor coupling f 2 R, the only coupling relevant in high- 
energy processes. (The couplings are defined later, in eqs. (J1J, (J5])). Early LHC data [22] 
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on measurements of W polarization in tt production and t-channel single-top production 
allow a limit of [—0.6, 0.3] to be put on f 2 R [IT] . There are stringent indirect constraints 
coming from low-energy measurements of anomalous tbW couplings. The measured rate of 
b — > sj puts stringent constraints on the couplings fm and f 2 L of about 4 x 10~ 3 , since their 
contributions to B- meson decay get an enhancement factor of mt /nib [231 El]. The bound 
on the anomalous coupling f 2R , however, is weak, viz., [—0.15, 0.57] at the 95 % CL [21]. In 
[25] a slightly more stringent bound has been found utilizing B& >s — B^ s mixing. It is thus 
clear that better direct limits on f 2 R are needed. 

We propose making use of the same lab-frame variables which we employed in [20] for 
tW~ production, but for a study of CPV, we would also need analogous variables for tW + 
production. The difference (in certain cases the sum) of these variables for the cases of 
tW~ and iW + production would be a measure of CPV. For convenience, we will adopt a 
linear approximation for the couplings. This enables determination of limits on the CP- 
violating parameter which are independent of the parameter itself. Moreover, we find that 
for maximal CPV, which we define in terms of equality in magnitude of and i couplings, 
there is no quadratic contribution to the asymmetry, making the linear approximation an 
excellent one. 

At the LHC, the initial state being pp is not self-conjugate and it is not evident that a 
study of CP-conjugate final-state particles can reveal the extent of CPV, if present. Never- 
theless, associated tW~ production is initiated by partons b and g, whose respective densities 
in p are equal to those of the conjugates b and g, from simple charge conjugation invariance 
of strong interactions governing parton distributions in hadrons. 

To see the consequences of CP invariance, consider a partial cross section for tW~ inclu- 
sive production at the LHC which may be written as 



da(p(p 1 )p(p 2 )-^t(p t , h t )W (pw)X)=Jdxidx 2 f b (xi)f g (x 2 )da bg ^ tw -(x 1 p 1 ,x 2 p 2 ,p t , h t ,p w ) (1) 

where da bg ^ tw - is the corresponding parton-level partial cross section, and f b , f g are the 
densities of b, g partons in the proton, pi, p 2 , p t and pw are respectively the momenta of the 
two protons, t and W~, and h t is the helicity of t. One can write an analogous expression 
for the partial cross section da (p{pi)p{p2) — > i(pt, hf)W + (pw)X > ) for iW + production. CP 
invariance at the parton level, which implies 

da bg ^ tw -(xip 1 ,x 2 P2,Pt,h t ,Pw) = da- bg ^ m +(x l pi,x 2 p 2 ,p u -h t ,p w ), (2) 

gives, for the hadron-level cross sections, 

da(pp -> t(p t , h t )W-{p w )X) = da(pp -> i(p t , -h t )W + (p w )X). (3) 

A violation of this relation would signal CPV. Thus, CP invariance implies equal and opposite 
longitudinal polarizations for t and t in the two processes. 

2. Effective vertices: We can write the most general effective vertices up to mass dimension 
5 for the four distinct processes of decay, t decay, production and i production as 

V t ^ bW+ = -(g/V2)V tb [Y(flLPL + flRPR)-^(cT^qJm w )(f 2L P L + f 2R P R )}, (4) 
V^ bW - = -( 9 /V2)V; b [Y(fiLPL + fmPR)-t(cr^q„/m w )(h L P L (5) 

Vb*-+tw- = -(0/V2)^[7^(/iVi + ^^ (6) 

v b ^ m+ = -G?/V2)^ 7 n/^+/r^ m 
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where f 1L , f 1R , / 2i , fiR, fiL, fin, /k, and / 2JJ are form factors, P L , P R are left-chiral and 
right-chiral projection matrices, and q represents the W + or W~ momentum, as applicable 
in each case. In the SM, at tree level, fn = fn = 1, and all other form factors vanish. New 
physics effects would result in deviations of fn and fn from unity, and nonzero values for 
other form factors. 

At tree level, i.e., in the absence of any absorptive parts in the relevant amplitudes giving 
rise to the form factors, the following relations would be obeyed: 

JlL — J1L, JiR — Jl-R) J2L — J2R, J2R ~ J2L- \p) 

This can be seen to follow from the fact that in the absence of final-state interactions 
leading to absorptive parts, the effective Lagrangian is hermitian. Thus, the hermiticity of 
the Lagrangian 



Cj = Hg/^)[VtJ>{l»UlLPL + hRPR)-(^ U /m w )d u W;U2LPL + f2RPR)}t 

+V£t{f (f 1L P L + JxrPr) - {a^/m w )d v W+ (f 2L P L + f 2R P R ) } b] + H.c. 



(9) 



from which the and t decay and production amplitudes d3J-(]7D m ay be derived, implies the 
relations flS}. 

On the other hand, if CP is conserved, the relations obeyed by the form factors are 

flL = flL', flR = flR', f2L = $2R] f2R = fiL- (10) 

We can write the phases of the form factors as sums and differences of two phases, one 
corresponding to a nonzero absorptive part, and another corresponding to CP nonconserva- 
tion: 

fn,R = \fiL,R\ exp (ia 1L:R + i5i L>R ); fi L;R = \ fi L;R \ exp {ia 1L , R - iSi L>R ); 

f2L,R = \f2L,R\ exp (ia 2 L,R + ^2L,r)] f%R,L = \f2L,R\ exp (iot 2 L,R ~ ^2L,r), 

where, for convenience, magnitudes of the couplings are assumed to be equal in appropriate 
pairs. In case there are no absorptive parts (a, = 0), we get the relations (JS]) as a special 
case. In case of CP conservation, we get the relations (fTUj) as a special case. In what follows, 
we will deal with a more general case when both CPV and absorptive parts are present. 
When calculating differences of variables for tW~ and iW + production, we will encounter 
only the combination 

A 2R = (Re/ 2 L - Re/ 2jR ) = -2\f 2R \ sin a 2 R sin 5 2 r. (12) 

(We assume fu = fn = 1 and V t b = 1 throughout). For A 2R to be nonzero, both absorptive 
parts and CPV have to be present, as can be seen from eq. (1121) . 

The effective vertices (HI)-© contribute to the production and decay of the top quark 
and antiquark. The parton-level Feynman diagram for tW~ is shown in Fig. (TJ with the 
effective vertex shown as a filled circle. The effect of anomalous couplings would show up 
in the production process in the total cross section, t or i angular distributions and or 
i polarizations. This is apart from W polarization, which we do not consider here. As for 
the decay process, for a given production angle and polarization of t/t , anomalous couplings 
would change the angular distribution of the decay products. 

In what follows, we obtain the contribution of CP-violating tbW couplings in sums or 
differences of observables in tW~ and iW + production and decay, and examine how well the 
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Figure V. Feynman diagrams contributing to associated tW production at the LHC The 
effective tbW vertex is shown filled circle. 

couplings can be constrained by them. We retain only Re/2_R and Ref 2 L , since these are 
the only ones which contribute in the limit of vanishing bottom mass, which we set to zero. 
We will assume that the anomalous couplings are small and a linear approximation of the 
couplings found to be good for |Re/2i?| up to about 0.05 [20J. 

3. CP-violating asymmetries: We now formulate CP-violating asymmetries in the single-top 
production subprocess gb — > tW~ and its CP-conjugated subprocess gb — > tW + through 
CP-violating anomalous tbW couplings. In [20] , we have obtained analytical expressions for 
the spin density matrix for the top quark produced in the process gb — > tW~ and for the 
top decay including the full contribution of anomalous tbW couplings. The corresponding 
expressions for the CP-conjugated process gb —> tW + and t — > bW~ can be obtained from 
those expressions by replacing Re/2i? by Ref 2 i taking appropriate helicity of the antitop. 

The simplest asymmetry we consider is the rate asymmetry in the production of tW~ 
and tW + . Assuming the equality of b and b densities in the proton, we can write, at linear 
order in the anomalous couplings, 

a(pp^tW-X) = a + Ref 2R a u a{pp -> tW + X) = a + Re/ 2 L cr x , (13) 

where <Jq is the SM cross section, which is identical for pp — > tW~X as well as pp — > tW + X 
and o~\ is the cross section arising from interference of the anomalous contribution with 
the SM amplitude, for unit anomalous coupling. We can then write the fractional rate 
asymmetry, 

a cpv = i/(2 < 7 )(Re/ 2i -Re/ M ) a x = l/(2a )A 2R a x . (14) 
In the linear approximation, the t polarization P t and the i polarization Pi can be written 

as 

P t = P t ° + Ref 2R Pt- P- t = -P t ° - Ref 2L P t \ (15) 

where P t ° and P] are contributions from the SM and from the pure anomalous couplings 
at linear order, respectively. In the linear order in anomalous couplings, a CP-violating 
asymmetry A^ v in top and antitop polarization defined by A^ v = Pf + Pt takes the form 

A^ v = -(Ref 2L - Ref 2R )Pt = -A 2K V\. (16) 

The asymmetries so far did not take into account top decay. A measurement of these 
asymmetries would need full reconstruction of the top. Thus the corresponding efficiencies 
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would have to be taken into account. We now look at asymmetries where an explicit decay 
channel of or t is taken into account. 

A rate asymmetry may be found corresponding to the combined production and semilep- 
tonic decay of the t and i, where both the production and decay processes get contributions 
from the anomalous couplings. Using the narrow-width approximation, we can write the 
cross section for the process pp — > tW~X — > bl + vW~ as 

a i+w - = a(pp bl + vW~X) = (l/r t ) x \p(pp -> tW~X) g> T(t ->• M+u)] , (17) 

where T t is total decay width of the top, p(pp —> tW~X), T(t — > bl + v) are respectively 
the production and decay density matrices (with appropriate integration of the phase space 
carried out), and (g> denotes a matrix product. The density matrix formalism is used here to 
ensure proper spin coherence between production and decay. Writing an analogous expression 
for the cross section <j£-w+ f° r the conjugate process, we define the asymmetry in charged- 
lepton production rates as 

Af PV = [cr i+w - - (j £ - w +]/[cr i+w - + a e - w +}. (18) 

Af gets contribution not only from the asymmetry in jt production, but also from the 
anomalous coupling in the decay into the leptonic channel. 

In all the following cases, we write the generic asymmetry in the tW~ process as as A~ , 
and that in iW + production as A + . We can then write, in the linear approximation, 

A- =A° + Ref 2R A 1 ] A + = A° + Ref 2L A\ (19) 

where A is the asymmetry in the SM, and A 1 is the contribution coming from the interference 
between the term with anomalous couplings and the SM term, for unit value of the coupling. 
The CP- violating asymmetry is then 

a cpv = A +_ A - = Am A \ (20) 

The azimuthal distributions of the charged lepton and the b quark arising from the top 
decay are sensitive to anomalous couplings [20] and we define the corresponding azimuthal 
asymmetries as 

A 4>e, h = W( C0S( Pi,b > 0) - ct(cos0£ i6 < O)]/[a(cos0£ ife > 0) + a (cos fa tb < 0)], (21) 

where fa, fa are respectively the azimuthal angles of the charged lepton and the b quark 
w.r.t. the top production plane. An analogous asymmetry can be defined for the charged 
lepton and b arising in antitop production. The CP-violating asymmetries A^ PV and 
are then defined as the differences between the relevant azimuthal asymmetries from the 
tW~ process and the conjugate process. 

The charged- lepton energy distribution is also sensitive to anomalous tbW couplings [20J, 
and the distributions for SM and different anomalous couplings all intersect at about Ef = 62 
GeV. To quantify the differences in distributions for different couplings, we constructed in 
[20J an asymmetry around the intersection energy Ef of the curves, defined by 

A Ee = [a(E e < Ef) - a(E e > Ef)\/[a(E l < Ef) + a(E e > E?)], (22) 

Using this, and the corresponding asymmetry for antitop production and decay, we define a 
CP-violating asymmetry Ag PV as their difference. 



6 



0.03 
0.02 
> 0.01 

Q_ 

° _ 

LU 

< -0.01 
-0.02 
-0.03 

-0.04-0.02 0.02 0.04 

A 2R 

Figure 2: The CP-violating energy asymmetry of charged leptons from t and i decay as a 
function of A 2 r for cm energy 14 TeV. The band corresponds to the la error interval in the 
SM for integrated luminosity 10 fb _1 . 



3. Numerical Results: We now evaluate the various asymmetries denned. For our numerical 
analysis, we consider for the LHC a cm energy y/s = 7 TeV for integrated luminosity L = 5 
fb -1 , and y/s = 14 TeV with L = 10 fb _1 . While the present run of the LHC is at y/s = 8 
TeV, we have checked that our results for this energy would be similar to those for y/s = 7 
TeV. We use the leading-order parton distribution function sets of CTEQ6L [26] with a 
factorization scale of mt = 172.6 GeV. Values of parameters used are: My/ = 80.403 GeV, 
a em {m z ) = 1/128 and sin 2 Q w = 0.23. 

We have imposed acceptance cuts on the charged-lepton rapidity and transverse momen- 
tum, viz., \r}\ < 2.5, and > 20 GeV. We have assumed an ideal situation of 100% efficiency 
for W and 6-jet identification, since we have used analytical expressions integrated over the 
full kinematic range for W and b. We later discuss the effect of realistic cuts needed for 
isolating the Wt final state. In order to discriminate the charge of the lepton which comes 
from top decay, we thus restrict ourselves to leptonic decay of the top and hadronic decay 
of the W. 

In the linear approximation, the asymmetry is given by the parameter A2R multiplied 
by the quantity A 1 for each asymmetry which is given in Table „ against the corresponding 
asymmetry. As can be seen, the lepton energy asymmetry is the largest among the asymme- 
tries constructed out of decay distributions. Also shown in Table 1, for various asymmetries, 
are the possible la limits on the CP-violating parameter A 2 r, using for the statistical un- 
certainty in A CPV the value 5A CPV = \/2/ y/LasM, where asm is the SM cross section for 
the relevant final state. For this, we have taken only one of the leptonic decay modes for the 
top, while including all hadronic channels for the decay of the W. The results are shown in 
Table 1. For the estimation of limits using A c / V and A C P PV , we assume 100% efficiency in 
the detection of the top, as well in the measurement of its polarization. The corresponding 
idealized limits are for comparison with the remaining more realistic ones obtained from 
kinematics of decay products. 

We see from the Table 1 that la limits on |A 2 rJ possible for y/s = 7 TeV with L = 5 fb _1 
are of the order of about 0.1, perhaps at the limit of validity of our linear approximation. The 
limit from A^> PV is nominally better, but cannot be realized since the polarization cannot 
be measured with 100% accuracy. For y/s = 14 TeV, on the other hand, the limits are of 
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A UFV for Am = 1 

Limit (L = 5 fb" 1 ) 


0.508 
0.0183 


0.526 
0.0180 


0.362 
0.106 


0.198 
0.175 


0.464 
0.075 


0.093 
0.127 


14 TeV 


A c^v for AaR = j 
Limit (L = 10 fb" 1 ) 


0.502 
0.00542 


0.536 
0.00506 


0.369 
0.0298 


0.186 
0.0343 


0.438 
0.0228 


0.086 
0.0386 



Table 1: The value of various CP-violating asymmetries for unit value of the CP-violating 
combination of couplings, A 2 r, and the corresponding la limits on |A 2 r| possible at the 
LHC. 



the order of a few times 10 -2 , well within the validity range of the approximation. In Fig. [2] 
we show A% PV , which gives the best limit on A2R, as a function of A2R. Also shown in the 
figure is the la band corresponding to the statistical uncertainty. The limit on |A 2 r| which 
would come from the measurement of v is 2.28 x 10~ 2 . 

We now discuss how the results of Table 1 would be affected on inclusion of realistic cuts 
needed for the isolation of the tW events and suppression of background. For this we make 
use of the analysis carried out recently by the ATLAS and CMS collaborations for the LHC 
data. 

We determine the factor by which the theoretical number for the Wt events get reduced 
by cuts in a single leptonic channel in top as well as W decay, making use of Table 2 of [13] 
and also Table 2 of [H]. We find this factor to be approximately 0.1. These searches use 
leptonic decay channels of W whereas a search for CPV should necessarily use the two-jet 
decays of W. For such a final state, ref. [5] indicates a somewhat better efficiency than 0.1 
for each channel, so we are being somewhat conservative. 

We then assume that the same factor can be used in our CP-asymmetry analysis (also 
at 14 TeV) to get an idea of the number of events which survive the cuts. Undoing the 
leptonic cut already used in Table 1, and including all leptonic channels (we assume detection 
efficiency of 0.6 for r's), we arrive at a correction factor of 0.325 for the event rates, and 
1.754 for the limits in Table 1. 

We thus see that in the best scenario of -y/i = 14 TeV with L = 10 fb -1 , a realistic limit 
of about 0.04 on |A 2 r| should be possible using the leptonic energy asymmetry. 

4- Conclusions and Discussion: We have investigated the possibility of measuring CP- 
violating tbW couplings at the LHC through the conjugate processes of tW~ and iW + 
production. We proposed a number of CP-violating asymmetries which would be sensitive 
to the CP-odd combination of couplings A 2 r = ]il — / 2 r, the difference of the couplings 
associated with the top and the antitop. 

We conclude from our analysis that the tW mode of single-top production is a good 
alternative process to look for CP- violating tbW couplings apart from a comparison of t and 
i decays in top-pair production. We find that the energy asymmetry, A^ v , is the most 
sensitive to the CP- violating parameter and would enable a limit of about 0.04 to be placed 
on |A 2 r|, for a cm energy of 14 TeV and integrated luminosity of 10 fb -1 . The limits possible 
for an operational energy of 7 TeV for integrated luminosity of up to 5 fb -1 are not as good, 
being at the level of 0.1-0.2. 

Higher-order QCD and electroweak effects can give rise to a partonic-level forward- 
backward asymmetry at the per cent level in the background process of ti production. This 
asymmetry can, because of valence-sea PDF difference, induce in the lab frame observables 
a charge asymmetry which does not originate in CPV. However, this is present only in the 



8 



sub-dominant quark-antiquark process. Moreover, how our leptonic observables would be 
affected needs detailed investigation, beyond the scope of this work. Table 1 shows that 
A% PV is fairly close in magnitude to A^ v , and seems to be a good measure of the polariza- 
tion asymmetry. Since polarization can only arise from chiral couplings, higher-order QCD 
effects cannot generate it. A% PV , at least, is therefore unlikely to get significant contribution 
from higher-order QCD effects. Top polarization due to electroweak effects is tiny [27]. 

All in all we conclude that it would be possible to obtain limits on the CP-violating 
anomalous coupling in the region of 10 -2 by employing the energy asymmetry of the charged 
lepton. Since this is the first investigation in this direction, these numbers corresponding to 
direct limits should be taken as encouraging. It would be worthwhile carrying out a more 
detailed and refined analysis. 

Various extensions of the SM would have specific predictions for these anomalous cou- 
plings. For example, the contributions to these form factors in 2HDM, MSSM and top-color 
assisted Technicolor model have been evaluated in [2S]- The predictions in these models, 
especially for the imaginary parts of the couplings, seem to be much smaller than the limits 
we expect from LHC with L = 10 fb _1 . 

Acknowledgment: SDR acknowledges financial support from the Department of Science 
and Technology, India, under the J.C. Bose National Fellowship programme, grant no. 
SR/SB/JCB-42/2009. 
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